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Abstract 

Previous studies indicate that lumbar radicular pain following disc herniation may be associated with release of several pro- 
inflammatory mediators, including interleukin-1 (IL1). In the present study, we examined how genetic variability in ILIA 
(rsl 800587 OT), IL1B (rsl 143627 T>C) and IL1RN (rs2234677 G>A) influenced the clinical outcome the first year after disc 
herniation. Patients (n = 258) with lumbar radicular pain due to disc herniation were recruited from two hospitals in Norway. 
Pain and disability were measured by visual analogue scale (VAS) and Oswestry Disability Index (ODI) over a 12 month 
period. The result showed that patients with the ILIA T allele, in combination with the ILIRN A allele had more pain and a 
slower recovery than other patients (VAS p = 0.049, ODI p = 0.059 rmANOVA; VAS p = 0.003, ODI p = 0.050 one-way ANOVA at 
12 months). However, regarding the IL1B/IURN genotype, no clear effect on recovery was observed (VAS p = 0.175, ODI 
p = 0.055 rmANOVA; VAS p = 0.1 05, ODI p = 0.21 4 one-way ANOVA at 12 months). The data suggest that the ILIA T/IL1RN A 
genotype, but not the IL1B J/IL1RN A genotype, may increase the risk of a chronic outcome in patients following disc 
herniation. 
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Introduction 

Studies of patients with lumbar disc herniation verified by MRI 
suggest that mechanical compression of the nerve-roots may 
induce lumbar radicular pain. However, disc herniation is also 
associated with an inflammatory response, with increased levels of 
pro-inflammatory mediators, such as interleukins (ILs), matrix 
metalloproteinases (MMPs) and prostaglandins (PGs) [1,2]. These 
inflammatory agents may promote disc degeneration, as well as 
sensitize the primary afferent sensory fibers. 

Several genetic variants may affect such processes and influence 
pain sensitivity. These might involve genetic variants important for 
the expression of GTP cyclohydrolase close to the nerve roots [3] , 
MMP1 in the intervertebral disc [4] or genetic variants in genes 
relevant to immunological systemic responses such as HLA [5] . In 
addition, genetic variability in ILIA may reduce the pain 
threshold after disc herniation [6]. 

The IL1 gene family includes ILla and ILlfi, as well as the 
interleukin-1 receptor antagonist (ILIRa). The former two are 
strong inducers of inflammation and one of the first cytokines 
produced by stress. They are secreted by a variety of cells, 
including monocytes, endothelial cells and disc cells [1,7]. The 
actions of these cytokines are, however, balanced by the expression 



of the endogenous ILIRa (encoded by IL1RN). This naturally 
occurring antagonist binds to IL1 receptors in competition with 
IL1, but does not elicit an intracellular response upon binding. 

ILla, ILlfi and ILIRa and their composite genotype have 
previously been linked to disc degeneration and low back pain in 
the general population [8,9]. We hypothesized that patients with a 
high expression of ILla (carriers of rsl 800587 T allele) or high 
expression of ILlfS (carriers of rsl 143627 T allele), combined with 
a reduced expression of ILIRa (carriers of rs2234677 A allele), 
would have an increased risk of long-lasting clinical pain as well. 
Hence, the purpose of our study was to investigate the combined 
effect ILIA, or IL1B, and IL1RN genotypes in recovery after disc 
herniation. 

Materials and Methods 

Patients 

Patients aged 18-60 years, with confirmed lumbar disc 
herniation by magnetic resonance imaging (MRI) with corre- 
sponding sciatic pain and positive Straight Leg Raising (SLR) test 
were recruited from Oslo University Hospital, Ullevaal, Norway 
and Haukeland University Hospital, Norway, for details see [10]. 
A total of 258 patients, all European-Caucasians referred to the 
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hospitals in 2007-2009, were included. However, at inclusion, 6 
patients changed their mind and did not want to participate, 
which gave data from 252 patients. In addition, 21 patients (8%) 
were lost during the follow-up. 

Ethics 

All participants received written information and signed an 
informed consent form. The study was approved by the 
Norwegian Regional Committee for Medical Research Ethics 
and the Norwegian Social Science Data Services. 

Clinical procedure 

After inclusion, the patients were followed up at 6 weeks, 6 
months and 12 months. Conservative treatment was received by 
42%, while the remaining 58% received surgery. At inclusion all 
patients underwent a standardized clinical examination with 
assessment of sensory and motor function and tendon reflexes at 
the lower limbs as well as an MRI scan. At 6 weeks and 1 2 months 
follow-up the clinical examination was repeated, while at 6 months 
follow-up patients reported their condition by a telephone 
interview and answered questionnaires by mail. 

All patients were asked to rate their pain intensity in activity 
during the last week on a 1 0-cm visual analogue scale (V AS) with 
Midpoints "no pain" and "worst possible pain". The validated 
Norwegian version of the Oswestry Disability Index (ODI) was 
used to assess problems with physical function related to low back 
pain [1 1]. The sampling of the clinical data was completed before 
the genotyping of the patients was performed. 

Genotyping 

Genomic DNA was extracted from whole blood cells using 
FlexiGene DNA isolation kit (Qiagen, Hilden, Germany). SNP 
genotyping was carried out using predesigned TaqMan SNP 
genotyping assays (Applied Biosystems) tor ILIA rs 1800587 OT, 
IL1B rsl 143627 T>C and IL1RN rs2234677 G>A. Approxi- 
mately 10 ng genomic DNA was amplified in a 5 (J.1 reaction 
mixture in a 384-well plate containing lx universal TaqMan 
master mix and lx assay mix, the latter containing the respective 
primers and MGB-probes. The probes were labeled with the 
reporter dyes FAM or VIC at 5 'end to distinguish between the two 
alleles. The reactions were performed on an ABI 7900HT 
sequence detection system (Applied biosystems) at the following 
program: After initial denaturation and enzyme activation at 95°C 
for 10 min, the reaction mixture was subjected to 60 at 95°C for 
15 s and 60°C for 1 min. Negative controls containing water 
instead of DNA were included in every run. Genotypes were 
determined using the SDS 2.2 software (Applied Biosystems). 
Approximately 10% of the samples were re-genotyped and the 
concordance rate was 100%. 

Statistics 

VAS activity score and ODI measurements over time were 
compared regarding genotypes by repeated measures analysis of 
variance (rmANOVA). When sphericity assumption was not met, 
a Greenhouse-Geisser correction was applied. Separate analyses 
were performed to check for potential confounding effects of the 
covariates age, gender, treatment and smoking status. Covariates 
with a p-value less than 0.1 were kept in the final model (Table 1). 
Further, VAS activity score and ODI score at 1 2 months after disc 
herniation were examined regarding genotypes by one-way 
ANOVA and Tuckey honestly significance difference (HSD) 
post-hoc comparisons. 
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Figure 1. The time course of clinical outcome measures following disc herniation. A) and B) patients grouped by ILIA CAT and IL1RN G/A 
genotypes. VAS activity score (p = 0.049 rmANOVA; p = 0.003 one-way ANOVA at 12 months), ODI score (p = 0.059 rmANOVA; p = 0.050 one-way 
ANOVA at 12 months). C) and D) patients grouped by IL1B T/C and IL1RNG/A genotypes. VAS activity score (p = 0.1 75 rm ANOVA; p = 0.105 one-way 
ANOVA at 12 months), ODI score (p = 0.055 rmANOVA; p = 0.214 one-way ANOVA at 12 months). Data are shown as means ± SEM. 
doi:1 0.1 371 /journal.pone.01 07301 .g001 



The combined effect of ILIA or IL1B, and IL1RN genotypes 
were analyzed. ILIA C/T and IL1RN G/A genotypes were 
grouped in the following variables: 1) ILIA C/C and IL1RN G/ 
G, 2) ILIA */T or IL1RN */C, and 3) ILIA */T and IL1RN */ 
C. IL1B T/C and IL1RN G/A genotypes were grouped in the 
following variables: 1) ILlp 1 C/C and IL1RN G/G, 2) IL1B */T 
or IL1RN */C, and 3) IL1B */T and IL1RN */C. Statistical 
analyses were performed by using the SPSS (version 20) statistical 
package. A p-value less than 0.05 was set as the level of statistical 
significance. The data are shown as mean ± SEM. 

Results 

The analysis of clinical outcome over time, i.e. from inclusion to 
12 months, revealed that the progression of pain and disability 
after lumbar disc herniation may be associated with the ILIA/ 
IL1RN genotype. In accordance with our hypothesis, patients 
with both rare alleles (ILIA */T and IL1RN */A) seemed to 
report more pain and have a slower recovery than other patients 
(Figure 1 A, B). A significant association between genotype and 
VAS activity score, and a borderline significant association 
regarding ODI score, were observed (VAS score jt> = 0.049, ODI 
score p = 0.059 within-subjects effects rm ANOVA). No clear 
association were, however, found for the clinical outcomes, with 
regard to the IL1B/IL1RN genotype (Figure 1 C, D; VAS score 
jf) = 0.175, ODI score p = 0.055 within-subjects effects rm AN- 
OVA). The characteristics of the cohort of patients stratified by the 
IL1A/IL1RN genotype and the IL1B/IL1RN genotype are listed 
in Table 2. 



Further analyses of the data sampled at 1 2 months showed that 
patients with the ILIA */T and IL1RN */A genotype had an 
increased one-year risk of a negative outcome, with more pain and 
poorer recovery (VAS score p = 0.003, ODI p = 0.050, One-way 
ANOVA). Post-hoc comparisons confirmed the differences 
between carriers of both ILIA T and IL1RN A alleles compared 
to none-carriers at 12 months (VAS score p = 0.006, ODI score 
p — 0.051, Tuckey HSD post-hoc comparisons). In contrast, no 
significant influence of the IL1B/IL1RN genotype on the one-year 
outcome was observed (VAS scored = 0.105, ODIp = 0.214, One- 
way ANOVA). Mean ± SEM values at 12 months are listed in 
Table 3. 

Discussion 

In the present study we investigated the clinical impact of the 
genetic variants ILIA rsl800587 C>T, IL1B rsl 143627 T>C 
and IL1RN rs2234677 G>A in patients following disc herniation. 
As expected, all patients suffered severe pain when they first were 
referred to the hospitals. However, patients carrying the ILIA T 
allele in combination with the IL1RN A allele (both rare alleles of 
these SNPs) reported more pain and disability during the follow- 
up period than the patients carrying only one or neither of these 
alleles. In line with previous data obtained from 131 occupation- 
ally active Finish men [8], our findings suggest that radicular low 
back pain may be dependent upon the ILIA T/IL1RN A 
genotype. 

The ILIA T allele has been associated with an enhanced 
promoter activity resulting in increased gene expression, both at 
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Table 2. Characteristics of patients grouped by the IL1A/IL1B and IL1RN genotypes. 







Gender, 

men/women (%) 


Mean age 
(min/max) 


Current smoker, 
yes/no (%) 


Treatment, 

surgery/conservative (%) 


ILIA C/T and ILIRN Gl A 


ILIA QIC and 
IL1RN G/G, n = 45 


20/25 (44/56) 


42 (22-59) 


12/33 (27/73) 


27/18 (60/40) 


ILIA 11* or IL1RN 
A/*, n = 167 


91/76 (54/46) 


40 (18-60) 


61/106 (37/63) 


98/69 (59/41) 


ILIA 11* and 
IL1RN A/*, n = 40 


24/16 (60/40) 


43 (25-59) 


19/21 (48/52) 


21/19 (52/48) 


ILJBJIC and ILIRN Gl A 


H.7B C/C and 
IL1RN G/G, n = 9 


4/5 (44/56) 


44 (32/60) 


3/6 (33/67) 


2/7 (22/78) 


/i!fiT/» or IL1RN 
A/*, n = 152 


85/67 (56/44) 


41 (18/60) 


54/98 (36/64) 


96/56 (63/37) 


/ZJBT/Cand 
IL1RN A/*, n = 91 


46/45 (51/49) 


41 (19/59) 


35/56 (38/62) 


48/43 (53/47) 



Min, minimum; max, maximum. 

doi:1 0.1 371/joumal.pone.OI 07301 .t002 



mRNA and at protein level, compared to the C allele [12]. 
Moreover, earlier data suggest that the IL1RN A allele, since it is 
in linkage disequilibrium with a VNTR polymorphism in intron 2 
of the IL1RN gene [13,14], is associated with reduced in vitro 
monocyte ILIRa synthesis and enhanced ILlfS production 
[15,16]. Thus, it seems likely that patients with the ILIA T/ 
IL1RN A genotype have overall enhanced levels of IL1 relative to 
ILIRa, resulting in an increased activation of IL1 receptors and a 
more pronounced inflammatory response. 

Analyses of intervertebral disc samples from patients suggest 
that a high IL1 to ILIRa ratio may be associated with 
degenerative disc disease [17]. IL1 may induce up-regulation of 
matrix-degrading enzymes, such as MMPs, and inhibit resynthesis 
of proteoglycans [17-19], thereby promoting degradation of the 
intervertebral disc matrix. Moreover, ILIA together with MMP3 
gene variation has previously been associated with Modic changes 
[20]. Genetic variability in MMP1 has also been suggested to 
contribute to low back pain and sciatica [4]. 

IL 1 , along with IL6 and TNF, secreted by a herniated disc may 
sensitize peripheral nociceptors. For example, ILlot increases the 
production of PGE2 in human disc cells in a dose dependent 
manner [2] , and alter the expression of substance P in rat DRG 



Table 3. Pain- and disability ratings at 12 months. 







VAS 


ODI 


ILIA C/T and ILIRN Gl A 


ILIA C/C and ILIRN G/G 


1.97 ±0.36 


13.10±2.08 


ILIA 11* or ILIRN M* 


2.30±0.20 


13.68±1.05 


ILIA 11* and ILIRN M* 


3.81 ±0.56 


19.68±2.86 


ILIB TIC and ILIRN Gl A 


\L1B C/C and ILIRN G/G 


2.13±0.93 


14.75 ±4.63 


/iIST/*or/iIRNA/* 


2.20±0.21 


13.25±1.17 


ILIB 11* and ILIRN A/* 


2.97 ±0.31 


16.63±1.52 



The table shows the 12 months VAS and ODI scores for the patients grouped by the combinations of ILIA C/T, ILIB T/C and ILIRN G/A genotypes. Mean ± SEM values 
are shown. 

doi:1 0.1 371 /joumal.pone.01 07301 .t003 



neurons [21]. ILlot is also active as a precursor. Peripheral and 
central administration of ILIRa has been found to alleviate 
inflammatory hyperalgesia in mice [22]. Hence, IL1 may 
indirectly facilitate nociceptive transmission, which may explain 
the increased pain experience reported by patients with the ILIA 
T/IL1RN A genotype. 

Regarding the ILIB T/ C SNP, it has been suggested that a shift 
from T to C results in a disruption of the TATA box resulting in a 
less active promoter. The T allele confers higher expression of the 
ILip gene compared to the C allele [23]. Actually, ILIB 
polymorphisms have previously been associated with symptomatic 
lumbar disc herniation [24], and together with ILIRN polymor- 
phism associated with low back pain in the general population [8] . 
In the present study, however, we did not find any associations 
between clinical pain and the ILIB T/ ILIRN A genotype. 

In addition to genetic influences, other risk factors such as 
physical work factors, psychosocial aspects and smoking may 
contribute to development of persistent back pain and disability 
[25-27]. Hence, the link between the IL1 genotype and persistent 
pain may be related to an increased inflammatory response, but 
also to more complex mechanisms. For example, the effect of 
genotype may be influenced by environmental factors [28]. 
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Moreover, a gene - lifestyle interaction is possible [29]. Our data 
suggested a possible association between smoking and persistent 
pain. However, both treatment and smoking were adjusted for in 
the final statistical analyses. 

In summary, our findings suggest that the development of 
chronic lumbar radicular pain may be associated with a high IL1 
activity relative to ILlRa. The combination oilLlA rsl800587 T 
allele and IL1RN rs2234677 A allele may increase the risk of a 
chronic outcome and persistent pain following lumbar disc 
herniation. However, the combination of IL1B rs 1143627 T 
allele and IL1RN rs2234677 A allele was not associated with 
increased pain in these patients. 
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